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We perform Monte Carlo calculation of correlation functions in 4d jV = 4 super Yang-Mills 
theory on R x S 3 in the planar limit. In order to circumvent the well-known problem of lattice 
SUSY, we adopt the idea of a novel large-rV reduction, which reduces the calculation to that of 
corresponding correlation functions in the plane-wave matrix model or the BMN matrix model. 
This model is a Id gauge theory with 16 supersymmetries, which can be simulated in a manner 
similar to the recent studies of the DO-brane system. We study two-point and three-point functions 
of chiral primary operators at various coupling constant, and find that they agree with the free 
theory results up to overall constant factors. The ratio of the overall factors for two-point and 
three-point functions agrees with the prediction of the AdS/CFT correspondence. 
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1. Introduction 

The gauge-gravity duality [jl|] has been one of the most important subjects in string theory over 
the past decade. The most typical example is the so-called AdS/CFT correspondence between type 
IIB superstring theory on AdS^ x S 5 and 4d JV = 4 U(A0 super Yang-Mills theory (SYM). Even in 
this case, however, a complete proof of the duality is still missing. In particular, the region described 
by the classical supergravity on the string theory side corresponds to the strongly coupled region in 
the planar large-Af limit on the SYM side. In order to study the strongly coupled 4d JV = 4 SYM 
from first principles, one needs to have a non-perturbative formulation such as the lattice QCD. 
The problem here is that the lattice regularization necessarily breaks translational symmetry, which 
is included in the supersymmetry (SUSY). In order to restore SUSY in the continuum limit, one 
generally has to fine-tune parameters in the lattice action. In fact any lattice formulations of 4d 
JV = 4 SYM proposed so far seem to require fine-tuning of at least three parameters 

Since 4d JV = 4 SYM has conformal symmetry, the theory on R 4 is equivalent to the theory on 
Rx S 3 through conformal mapping. The novel large-A^ reduction [^] connects the planar large-A^ 
limit of this theory to a reduced model, which can be obtained by shrinking the S 3 to a point. The 
resulting one-dimensional gauge theory with 16 supercharges can be studied by using the Fourier- 
mode simulation [Q] as in recent studies of the DO-brane system []|]. Thus we can perform Monte 
Carlo calculations in 4d JV = 4 SYM respecting SUSY maximally and without fine-tuning. 1 

In this article we present explicit results for correlation functions of chiral primary operators 
(CPOs) in 4d JV = 4 SYM. 2 In particular, we find that the two-point and three-point functions agree 
with the free theory results up to overall constant factors even at fairly strong coupling. Moreover 
the ratio of the overall factors agrees with the prediction of the AdS/CFT correspondence. 3 

2. Large-N reduction for JV = 4 SYM on R x S 3 

Let us first discuss the novel large-Af reduction for JV = 4 SYM on R x S 3 . By collapsing the 
S 3 to a point, we obtain the plane wave matrix model (PWMM) or the BMN matrix model [^] 4 , 
whose action is given by 



Spw = ~y~ I <ltlr 

<?pw 



l -{D t X M f - l -[X M ,X N ] 2 + Vd/F- VrfM,^] 



2. l 1 



(2.1) 



Here the parameter jx is related to the radius of S 3 as R S 3 = jj, and the covariant derivative is defined 
by D t = d t — /[A, • ], where A(t), as well as X M (t) and *P(f)> is an x hermitian matrix. The 
range of indices is given by 1 < M,N < 9, 1 < z, j, k < 3 and 4 < a < 9. The model has the SU(2|4) 
symmetry with 16 supercharges. 



'See refs. [n] for proposals for finite N. 

2 See ref. ||7|| for some preliminary results on the Wilson loop. 

3 There are also Monte Carlo studies of the 4d .,-¥ = 4 SYM based on matrix quantum mechanics of 6 bosonic 
commuting matrices [Q], which give results consistent with the AdS/CFT for the three-point functions of CPOs. 

4 Properties of this model at finite temperature are studied at weak coupling Jlo| , |T"i] ] and at strong coupling [ p^ j. 
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In fact the model possesses many discrete vacua representing multi fuzzy spheres, which are 
given explicitly by 

V v 

X i = n($)(L < f l ' ) ®l kl ) with Y,n I k I =N, (2.2) 

7=1 1=1 

(r) (r) (r) 

where L\ are the r-dimensional irreducible representation of the SU(2) algebra [L\ ,Lj ] = 

(r) 

i£ijkL k . These vacua preserve the SU(2|4) symmetry, and are all degenerate. 

In order to retrieve the planar jV = 4 SYM on R x S 3 , one has to pick up a particular back- 



ground from (2.2), and consider the theory (2.1) around it. Let us consider the case 



2 

and take the large-TV limit in such a way that 



n I = n + I---p- for /=!,••• ,V , (2.3) 



OO — 



" ^ -o , v -»• oo , with A PW = fixed . (2.4) 



v 

Then the resulting theory is claimed ^ to be equivalent 5 to the planar limit of jY = 4 SYM on 
Rx S 3 with the 't Hooft coupling constant given by 

^sym = 2tc 2 Xpw(Rs 3 ) 3 = ■ (2-5) 

n pL A 



The above equivalence may be viewed as an extension of the large-TV reduction p5Q, which 
asserts that the large-TV gauge theories can be studied by dimensionally reduced models. The orig- 
inal idea for theories compactified on a torus can fail due to the instability of the U(1) D symmetric 



vacuum of the reduced model [16]. This problem is avoided in the novel proposal since the PWMM 
is a massive theory and the vacuum preserves the maximal SUSY. Since the planar limit is taken 
in the reduced model, the instanton transition to other vacua and the "fuzziness" of the spheres 
are suppressed. Viewed as a regularization of the JV = 4 SYM on Rx S 3 , the present formula- 
tion respects the maximal SU(2|4) symmetry (with 16 supercharges) of the PWMM, and in the 
limit ( |2.4[ ) the symmetry is expected to enhance to the full superconformal SU(2,2|4) symmetry 
(with 32 supercharges). Since any kind of UV regularization breaks the conformal symmetry, this 
regularization is optimal from the viewpoint of preserving SUSY. 



3. Correlation functions in 4d J/ = 4 SYM 

As simple examples of 1/2 BPS operators, let us consider the CPOs given by 

<(*) = T ai ... aA tr (x*X#:.X*(xj) , (3.1) 

where T av .. aA is a symmetric traceless tensor and X* represents the six scalars in 4d Jf = 4 SYM 
on R 4 . Thanks to the conformal symmetry, the forms of two-point and three-point functions of the 

5 See refs. JTift for earlier studies that led to this proposal. This equivalence was checked at finite temperature in the 
weak coupling regime wM. It has also been extended to general group manifolds and coset spaces KM. 
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CPOs are determined as 

/ <(x 1 )<(x 2 )<(x 3 )) = c AlA2A3 (^( Xl )^(x 2 )^(x 3 ) 



I free 



free 



(3.2) 



where c A and ca!A 2 a 3 are over-all constants depending on Asym in general, and (• • • )f ree denotes the 
results of free theory. The analysis on the gravity side suggests [ |l7| ] 



C A[A2A3 



y ' A \C- ' A jt 'A3 



CA1A2A3 



N- 



yj C A \ C Al c A ^ 



1 for V A, . 



(3-3) 



In order to relate the above operators to those in the PWMM, we first perform the conformal 
mapping 6 from R A to R x S 3 . Then the M-point functions of the CPO ff^ s3 on R x S 3 are related 
to those in PWMM as 



dQ, 



(i) 



2tt 2 



dQ. 



(M) 



2% 2 



fRxS 



(3.4) 

where we have defined 0™(t) = T av .. aA tr (x ai X a2 ■ ■ -X aA (f)) [§. 

We calculate the two-point functions ^trZ^r^trZ 1 " 2 ^)^, where Z = ^(X 4 + /X 5 ), and the 
three-point functions (tr (X4X 5 (fi)^ tr (x 5 X 6 (t 2 )^ tr (x^fe)) J. The CPOs we consider here 

3 /2 

have A = 2, and the AdS/CFT predicts C222 = c 2 , which we test by Monte Carlo calculations. 



4. Monte Carlo method 



In order to simulate the PWMM (|2.1[), we compactify the £-direction to a circle of circumfer- 
ence j3. Since we are interested in the properties at zero temperature, we impose periodic boundary 
conditions on both scalars X(?) and fermions *P a (f), which keep SUSY intact. In Fourier-mode 
simulation [|J], we first fix the gauge symmetry completely by choosing A(t) = ^diag(ai , • • • , a#) 
with —71 < a a < K, and then make a Fourier expansion X,(r) = £„ = _ A X ;! „e ,Q " tf (ft) = ^) and sim- 
ilarly for the fermions. The upper bound A on the Fourier modes plays the role of the UV cutoff. 
The original PWMM can be retrieved by just taking the limits j3 — > °° and - — > 00 since there are 
neither UV nor IR divergences. The model regularized by finite j3 and A can be simulated by the 
RHMC algorithm. This method has been applied extensively to the DO-brane system corresponding 
to fi = 0, and the results confirmed the gauge/gravity duality for various observables [^]. 7 

Since the parameter gp W in the action (11) can be scaled away by appropriate redefinition of 



fields and parameters, we take gp W N = 1 without loss of generality as in refs. [||]. In this convention 



one finds from eq. ( \1.5\ > that the small (large) jj. region in the PWMM corresponds to the strong 



(weak) coupling region in the 4d jV = 4 SYM. 

6 The metrics of /? 4 and R x S 3 are related as ds^ = dr 2 + r^dQ.^ = e^'ds^^ , where r=^e^ r . The transformation 

of the CPOs is given by = e^' &f . 

7 See refs. wM for Monte Carlo calculations based on the lattice regularization. 
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Figure 1: (Left) The two-point function (trZ 2 (p) \xZ^ 2 (—p)j is plotted in the log scale. The curves repre- 
sent the corresponding free theory results multiplied by 0.919, 0.799, 0.647 for/x =4.0,2.0, 1.3, respectively. 
(Right) The ratio of the two-point function to the corresponding free theory result is plotted in the linear scale 
for fi =4.0,2.0,1.3. 




Figure 2: (Left) The three-point function ^tr (x 4 X 5 (p)^ tr (x 5 X 6 (0)j tr \ X 6 X 4 (-p)j \ is plotted in the 
log scale. The curves represent the corresponding free theory results multiplied by 0.850, 0.716, 0.491 for 
ji = 4.0,2.0, 1.3, respectively. (Right) The ratio of the three-point function to the corresponding free theory 
result is plotted in the linear scale for jj, — 4.0,2.0, 1.3. 



5. Numerical results 

The parameters describing the background ( |2.3| ) are chosen asn=|,V = 2,£ = 2, which 



corresponds to the matrix size N = 6. We use ( |2.2[ ) with ( |2.3| ) as the initial configuration and check 



that no transition to other vacua occurs during the simulation. The values of }X we use are /I = 
4.0,2.0, 1.3, which correspond to Asym — 0.55,4.39, 16.0, respectively, in the chosen background. 
Thus we cover a wide range of the coupling constant. The regularization parameters in the t- 
direction are taken as j3 = 5.0, A = 12 for all cases. 

In fig. U (Left) we plot the two-point function 8 (trZ 2 (p)trZ^ 2 (—p)^. We find that the results 
agree well — up to overall constants depending on jj. — with the corresponding free theory results, 



s The Fourier transform of an operator Git) is defined as S'(p) = j [q dt ff(t) 
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which are obtained by just switching off the interaction terms in the reduced model with the 
regularization parameters. In fig. [j] (Right) we plot the ratio to the free theory results. Fig 
shows similar results for the three-point function (tr (x^X^{p)^ tr ^X 5 Xg(0)^ tr (X(X^{— p) 

We can extract the the overall constants cor- 
responding to C2 and C222 in eq. (3.2) from figs, [j] 
and ||, respectively. Since the data on the right 
panels are not completely constant in p, we take 
the maximum and minimum values as the upper 
and lower bounds of the estimates. In fig. || we 
plot the overall constants obtained in this way 
for three values of pi. The data points represent 
the mean value of the upper and lower bounds. 
We find that our results for various coupling con- 
stants lie on the straight line which represents the 
prediction C222 



same 
ure 




,3/2 



Figure 3: The overall constants corresponding to q 

from the AdS/CFT. Our and c 2 2i in eq. (fO) are plotted in the log-log scale. 

3/2 



results therefore suggest that the relation (3.3) 
holds also at intermediate coupling constants. 



The straight line presents the relation C222 
predicted by the AdS/CFT. 



6. Summary and discussions 

We have made the first attempt to investigate nonperturbative properties of the 4d JV = 4 S YM 
from first-principle calculations. Our formulation is considered optimal in preserving SUSY, which 
seems to give us the virtue of making the coupling constant dependence of the CPO correlation 
functions restricted mostly to the overall factors. This feature of our formulation enables us to test 



the prediction (3.3) of the AdS/CFT correspondence already for quite a small matrix size. Our 
results suggest that the relation extends to intermediate coupling constants as well. 

In fact there is strong evidence from field theoretical analysis in the 4d jV = 4 SYM that the 
non-renormalization theorem holds for two-point and three-point correlation functions of CPOs 
Jl9| , |20|] 9 , which implies ca = ca 1 a 2 a 3 = 1. It is therefore expected that the data points in fig. || will 



approach C2 = C222 = 1 as we take the limit (2.4) for fixed Asym> which needs to be checked. 



The analysis of four-point functions would be more interesting [E2] since there is evidence 



for the non-renormalization theorem only in the extremal and next-to-extremal cases [23], and in 
fact the AdS/CFT predicts its violation for the other general cases in the strong coupling limit J24|]. 
We consider it very interesting that such nonperturbative issues in 4d JV = 4 SYM have become 
accessible by computer simulations. 
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